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This is a study of a possible alternative procedure for adding a potential energy to the free electron Dirac 
equation. When Dirac added potentials to his free electron equation, there were two alternatives (here called Dl 
and D2). He chose Dl and lost charge conjugation symmetry, found Ehrenfest equations that depended on the 
sign of the energy of the state determining the expectation value, encountered Klein tunneling, zitterbewegung 
and the Klein paradox. The Dl alternative also predicted that deep potentials should pull positive energy states 
down into the negative energy continuum, possibly creating an unstable vacuum. Extensive experiments (1975- 
1997) found no evidence for this instability, but did find low energy electron-positron pairs with sharply defined 
energies and unusually low counting statistics. These pairs tended to disappear with higher beam currents. This 
paper explores the other alternative, here called D2 and finds charge conjugation symmetry preserved, Ehrenfest 
equations are classical, Klein tunneling is not present, unstable vacuua are forbidden, zitterbewegung is absent 
in the charge current density, new excitations of bound electron-positron pairs are possible in atoms, and the 
energies at which low energy electron-positron pair production in heavy ion scattering occurs is well described. 
Also all of the positive energy calculations, including those with the Coulomb potential, the hydrogen-like atom, 
are retained exactly the same as found in alternative D 1 . It might have been better if Dirac had chosen alternative 
D2. 



I. INTRODUCTION 

In the period 1973-1999 there was considerable activity 
seeking to use the attractive potential of two heavy ions scat- 
tering to pull down the lowest bound state electron states of 
the two nucleon system down through the mass gap (— m, to) 
, (here c = 1), into the negative continuum states. The ideafT- 
[3|] was to destabilize the Dirac Sea by bringing an electron 
state into contact with the negative energy states. The possi- 
bility of making the vacuum unstable generated a great deal 
of theoretical and experimental excitement which was docu- 
mented in three conference proceedings |]4J-|6|]. The electronic 
Dirac equation seemed to predict fascinating physics if the 
positive energy bound states could be pulled down into the 
negative continuum. 

However, after many experiments the only surprising 
observation||4j-|6|] was the existence of a very narrowly defined 
total energy peak for an electron-positron pair at a low energy 
of hundreds of KeVs. The best theory had predicted that 
there should be a dependence on the total charge Zi + Z2 
of the two heavy nuclei and that there should be a critical 
threshold for this sum. Above this critical sum there should 
be evidence of vacuum instability and nothing exceptional 
should be observed below this threshold. In fact, the electron- 
positron sum energies for various heavy ion pairs did not cor- 
relate well with the nuclear charges and sub-critical pairs |7] 
also displayed the sharply peaked total energy e^e+ pairs. 

The early experiments were capable of detecting only one 
member of the electron-positron pair, but as the experiments 
grew in sophistication, detection of both members of the pair 
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was possible and the narrowness of the total energy distribu- 
tion was verified. The mysterious source of these e^e+ was 
apparently generated in the collisions and moved at a slow 
speed from the scattering center The counts for these e^e+ 
pairs were generally not very large. They were difficult to 
detect and it took very long runs to generate good statistics. 

A collaboration called APEXfS] was organized through Ar- 
gonne National Laboratory to use the ATLAS [91 heavy ion 
accelerator, which was capable of large ion currents, and to 
use specially developed spectrometers for detecting simulta- 
neously the electrons and positrons from the collisions. The 
hope was to use the larger currents to lower run times for gen- 
erating adequate statistics. The early runs[10] at low beam 
currents generated fairly convincing data that a sharp total en- 
ergy peak was present, though the counting statistics were low. 

However, as the beam current was increased, the peak did 
not grow out of the background, rather it seemed to disappear 
as the run proceeded. The properly skeptical experimental- 
ists decided that the failure of the sharp peak to survive large 
beam currents and long running times indicated it was an ar- 
tifact, an unexplained, unreliable chimera lll lll . lll2ll .lll3ll. De- 
spite reasonable objectionsi 14|, the APEX results essentially 
closed down research in this whole area. There have been al- 
most no papers on this subject since 1999. 

What could have gone wrong? The Dirac equation for deep 
square wells predicts that bound state energy levels can be 
pulled below the bottom of the mass gap. Yet no vacuum in- 
stability was observed. None of a great variety of theoretical 
ideas could predict the low energy, sharply defined total en- 
ergy e^e+ peaks. Also none of the theories could explain ei- 
ther the wide spread observation of these e~e^ pairs for many 
heavy ion pairs at low beam currents or the disappearance of 
the sharp peak as the beam current increased. What could be 



the reason that the efforts of so many capable scientists were 
not able to confirm the predictions and explain the sharp e^e+ 
peaks? 

A very large number of theoretical ideas have been explored 
hoping to explain these failures. So far it appears that no one 
has gone back to the beginning of relativistic quantum me- 
chanics and considered the process that Dirac used to add in 
the potentials to the free electron theory. Before Dirac added 
in the potential energy, his free electron equation had charge 
conjugation invariance and the equations of motion of the ex- 
pectation values obeyed the Ehrenfest classical equations. 

The goal of this study is to re-examine Dirac's 
derivation lITsIl . III6I1 of his equations and the way he included a 
potential energy in his free electron equation. The possibility 
that there may have been a alternative "path not taken"|[17| 
is explored and this possibility seems to recover much of 
Dirac's results, resolves the above issues and also points to 
new physics that may have been excluded by Dirac's deriva- 
tion. In the following, the path taken by Dirac will be labeled 
Dl and the other alternative will be labeled D2. Before de- 
lineating these alternatives, it is useful to summarize some of 
the theoretical consequences of alternative Dl. 



A. Consequences of Dirac's Choice for Incorporating Potential 

Energy 

The non-interacting Dirac equation has an energy spectrum 
symmetric about zero ener gy a nd the equation obeys charge 
conjugation invariance lllSll . lll9ll . That is, the negative energy 
solutions are mapped onto the positive energy solutions and 
vice versa. Yet, when a 4-vector potential A^^ is included 
using alternative Dl, the charge conjugation invariance breaks 
down and the energy spectrum is no longer symmetric about 
zero energy. 

The sign of the vector potential is seen to change as the 
charge conjugation transformation is applied to solutions of 
the interacting Dirac equation. This has long been interpreted 
to mean that the positrons, which are represented by the trans- 
formed wave functions ^e have the opposite charge to that of 
the electrons flal, [181, fl91. However, the sign of the positron 
was originally assigned by looking at the change in the charge 
of the Dirac vacuum in the absence of a negative electron. 
This assignment of charge is valid even in the non-interacting 
Dirac equation where there is no vector potential or explicit 
charge in the Hamiltonian. This observation would support 
the argument that the sign change in the vector potential in 
the transformed Dirac equation could simply be interpreted as 
evidence that the charge conjugation invariance of the Dirac 
equation breaks down. The interpretation of the charge of a 
hole in the vacuum and its positive energy relative to the vac- 
uum is independent of the charge conjugation transformation 
itself. 

In the Dl alternative, Dirac's original approach, which so 
strongly informs our intuition, we have an apparently consis- 
tent theory of electrons and positrons, even though there are 
some inconsistencies: (1) failing to be charge conjugation in- 
variant in the presence of a vector potential, (2) having Ehren- 



fest equations of motion for expectation values that depend on 
the sign of the energy of the wave function forming the expec- 
tation value, and (3) having the existence of Klein tunneling 
in which low energy electrons appear to tunnel into high and 
wide potential steps that are higher than the mass gap of 2m. 
Finally, the perplexing failure of the heavy ion scattering pro- 
gram to discover an unstable vacuum and the appearance of 
narrowly defined total energy of e^e+ pairs, which seem to 
disappear at high beam currents, suggest that the Dl alterna- 
tive could be inconsistent with nature. 

The Ehrenfest theorem, that the equations of motion of ex- 
pectation values of various operators should reproduce the 
classical equations of motion 1. 1& ,19.] . has long been a lim- 
iting check on the validity of quantum mechanical theories. 
Derivations in the alternative Dl of the equation of motion for 
the canonical momentum for a particle (electrons with nega- 
tive charge) in an electric and magnetic field from the inter- 
acting Dirac equation ends up with the sign of the Lorentz 
force law depending on the sign of the energy of the state be- 
ing used to determine the expectation value of the canonical 
momentum lllSll . 

A similar Dl derivation from the interacting Dirac equation 
of the motion of a spin in a magnetic field shows that the di- 
rection of the torque on the spin again depends on the sign of 
the energy of the state being used to calculate the expectation 
of the spin componentsi 18|. Both of these results contradict 
the Ehrenfest theorem that the equations of motion for the ex- 
pectation values should be classical. 

Similarly, the probability current density in Dl should be 
a constant in the limit of no electric field, but it remains time 
dependent in the free particle limit while the momentum re- 
mains constant. This also represents a violation of the clas- 
sical limit of an Ehrenfest equation. The probability current 
density also exhibits zitterbewegung oscillations at high fre- 
quencies of 2m. 

The Klein paradox lll9ll . ll20ll has disclosed a further com- 
plication with the interacting Dirac equation in the Dl alter- 
native. There have been a wide variety of demonstrations 
that something strange occurs if a potential stepis too large in 
energy. Recently, Dragoman[21J and Bowen[22] have noted 
that the traditional view of the Klein paradox, that the reflec- 
tion and transmission coefficients are not positive, is easily 
resolved analytically, but the Klein tunneling remains unre- 
solved. 

Finally, combining these theoretical observations of the 
breakdown of charge conjugation invariance, defective Ehren- 
fest equations of motion, and Klein tunneling, with the ex- 
perimental failure to observe a straight forward prediction 
of vacuum instability and the lack of a theory for the sharp 
e~e+ pairs, argues that a second look at Dirac's procedure for 
adding in the potentials should be undertaken. 

The rest of the paper is organized as follows: In section 
II the two alternative approaches to relativistic quantum me- 
chanics will be delineated with the following: 

1. Review of Dirac's free particle equation and its solu- 
tions. 

2. Addition of the potential energy to the free electron 



equation and description of the Dl and D2 alternatives. 



A. Free Particle Equations 



Section III presents several topics comparing the alternatives 
Dl and D2 in the following sub-sections. 

1. Charge conjugation invariance survives the addition of 
vector potentials n D2, while it fails in Dl. 

2. Ehrenfest equations deficiencies in Dl and resolution 
of these in D2. 

3. Effect of different one electron states used in 131 and 
£12. 

4. Bound states in the mass gap in Dl and D2. 

5. Differences between Dl and D2 forpiecewise constant 
potentials. 

6. Klein tunneling in Dl and D2. 

7. Proof by contradiction of a limit for attractive potentials 
inD2. 

8. Aspects of Dirac hole theory in both alternatives. 

9. Probabililty current density and zitterbewegung in Dl. 

10. Charge current density and the absence of zitterbewe- 
gung in D2. 

11. Completeness in D2. 

12. Transformations and Feynman Stuckleberg theory in 
£11 and D2. 

13. Feynman diagrams and Perturbation Theory in D2. 

Section IV compares the consequences of D2 and Heavy 
Ion Experimental Data with the following subsections; 

1 . Historical Introduction 

2. Hydrogenic Bound states in the Mass Gap 

3. Decay modes for bound e^e+ metastable states 

4. Cross section as a function of energy 

5. Free e^e+ energies in the Ion Center of Mass Frame 

6. Free e^e+ energies transformed to the Lab Frame 

7. Comparison between Experimental Peaks and bound 
e~e+ metastable transitions 

8. Explanation of Tables 1 and 2 

Section V is a summary which is followed by two appen- 
dices deriving Ehrenfest equations and one discussing induced 
metastable state decays. 



n. DESCRIBING THE TWO ALTERNATIVES 

In this section the motivation for the D2 alternative and its 
distinction from Z31 is delineated. 



Dirac 1 1611 sought matrices a and /? so that the free particle 
Hamiltonian could be written (c ~ 1 and h — 1) for momen- 
tum p and mass m as 



Hq = a ■ p + I3m. 



(1) 



The eigenstates of this Hamiltonian were of two types: pos- 
itive energy u{p) and negative energy v{p) and each of these 
satisfied the equations 



{a ■ p + /3m)u{p) = Epu{p), 



{a ■ p + (3m)v{p) = —Epv{p), 



and where 



Ep = ^/jp+^. 



(2) 



(3) 



(4) 



One of the most important ideas that Dirac discovered was 
the charge conjugation symmetry that connected these two 
eigenvectors to each other and enabled an understanding of 
the positron. To briefly review this symmetry, one takes the 
negative energy eigenvalue equation and changes » — >^ — p 
and takes the complex conjugate of the equationplSI. This 
yields 

{-a* ■p + l3*m)v*{-p) = -EpV*{-p), (5) 

a 
matrix C that has the following properties 



where use of Ep = E^p has been made and we now seek a 



and 



Ca*C- 



C(3*C- 



(6) 



(7) 



Substituting these into Eqn. 5 yields 

{-cx-p-(3m)Cv*{-p) = -EpCv*{-p), (8) 

which becomes the positive energy eigenvalue equation af- 
ter factoring out —1. This means that within a phase factor 
we have the following connection between the two different 
energy-signed eigenstates 



uip)^Cv*i-p). 



(9) 



Thus, the energy spectrum is symmetric about zero. For ev- 
ery negative energy eigenvector v{p) at energy —Ep, there is 
a corresponding positive energy eigenvector at +Ep which is 
given by the charge conjugation transformation. In the next 
section the extension of these free particle results to include 
potential energies following Dirac alternative Dl and the al- 
ternative path D2 will be delineated. 



B. Incorporating Potential Energy 

Given the symmetries of the non-interacting Dirac equa- 
tion, the next question is how to add in a potential energy. 
Intuitively, for the positive energy eigenstates the inclusion of 
a small potential V should have the form 



E — En 



V, 



(10) 



and should give the correct non-relativistic limit as the mo- 
mentum goes to zero. 



E=^ 



V 



PL 

2m 



V. 



(11) 



For the negative energy eigenstates, there is no clear cri- 
terion for choosing the sign of the potential energy to be in- 
cluded. There are two possibilities that could make sense: 

Alternative Dl: (Dirac's choice) is simply to add the po- 
tential to the negative energy eigenvalue. 



E 



-Ep + V 



-y/p^+~^+V K -m- — + V. (12) 



2m 



Alternative D2: (The path not taken) would add the poten- 
tial with a negative sign 



E 



-Ep-V 



Vp^ + ' 



-Vx -m- — -V. (13) 
2m 



At this moment there is actually little to distinguish or jus- 
tify either alternative. The first alternative Dl is the sim- 
plest, applying the same rule to both energy signed eigenval- 
ues. The D2 alternative restores the symmetry about zero 
energy that the Dl alternative seems to destroy. 

The Dl alternative has already been shown to be rela- 
tivistically invariant under Lorentz transformations and leads 
to an immediate proof that the interacting Dirac equation is 
covariant lUQII . 

Is it possible that the D2 alternative can be made in a 
Lorentz invariant fashion that will allow a proof that the re- 
sulting Dirac equation is covariant? The first step in such a 
demonstration requires the review of the Casimir projection 
operators S+(p) and B^(v)(l^ 



and 



B+{p)u{p) = u{p), B+{p)v{p) = o (14) 



B- (pHp) = v{p), B_ {p)u{p) = 0. (15) 



The sign of the energy is a Lorentz invariant. The pos- 
itive and negative energy states do not mix under Lorentz 
transformations! 19]. Using these Casimir operators an op- 
erator that yields the Lorentz invariant, the sign of the energy, 
can be written as 



sgniEip)) = B+{p) - B^{p). 



(16) 



So, a Lorentz invariant way to write the D2 alternative 
is to substitute sgn{E{p))V for the potential energy in the 



Dirac equation. Then Dirac's equation with a potential en- 
ergy would have the form 



{a ■ p + (3m + sgn{E)V)tp{p) — E^{p), 



(17) 



which results in the positive eigenvalue equation being written 

as 

{a-p + (3m + V)u{p)=^ Eu{p), (18) 

and the negative energy eigenvalue equation being written as 

{CX-P + I3m-V)v{p) = -Ev{p). (19) 

If we again apply the charge conjugation transformation, 
we see that the energy spectrum in the presence of a potential 
continues to have the symmetry that was present for the free 
particle solutions. Charge conjugation transformation invari- 
ance is restored. 

If the more general case of coupling to electromagnetic 
fields is considered, the minimal coupling interaction for this 
D2 alternative in the Dirac equation should be replaced by 



pM -^ pt' - sgn{E)eA'' , 



(20) 



where the charge of the electron is e (e < 0) and we attempt 
to follow the sign conventions of Bjorkan and Drellfl9]. Be- 
cause the sgn{E) is a Lorentz invariant, there is no change in 
the vector nature of this transformation and the usual proof of 
covariance carries through automatically. 



III. COMPARISONS BETWEEN Dl AND D2 

In this section several properties of both the Dl and D2 
alternatives will be reviewed. Several rather striking differ- 
ences will be found and the reader is urged to keep in mind 
that this comparison will in some cases be jarring to our intu- 
itions, schooled as we are in the Dl alternative. 



A. Charge Conjugation Invarlance with potentials In D2 

The eigenvalue equation for the Dirac Hamiltonian in the 
D2 alternative will now have the form 

{a. ■ {p — sgn{E)eA) + /3to + sgn{E)e(j))il;{p) — Eipip)- 

(21) 
For a positive energy state we would have 

(a • (p - eA) + I3m + e(f))u{p) = Eu{p) (22) 

and for a negative energy state 

[a- [p + eA) + I3m~ e0)w(p) = -Ev{p), (23) 

where u{p) and v{p) are used to represent the positive and 
negative energy solutions in the presence of the vector poten- 
tial A^. Notice that the positive energy equation in D2 is ex- 
actly the same as in Dl, so, for example, the positive energy 



hydrogenic solutions with the Coulomb potential will be the 
same. It may be likely that a reader might confuse alternative 
D2 with an earlier alternative studied by Miiller, Rafelski, and 
Soff|23]. In this Miiller alternative the Coulomb potential was 
assumed to be multiplied by the j3 matrix. For a hydrogenic 
solution with a Coulomb potential this altered Dirac equa- 
tion can be solved analytically, but the atomic spectra does 
not agree with experiment. The alternative D2 is completely 
different from this earlier model. In D2 the positive energy 
hydrogen is completely unchanged and thus all of the atomic 
states that agree with experiment are recovered. 

If we now take the negative energy equation and change 
p — > —p and take the complex conjugate we obtain 

OL* ■ {-p + e A) + 13*171 -e(j))v*{-p) = -Ev*{-p). (24) 

Looking for a matrix C exactly as before, we obtain 

{-a- {p~ eA)- Pm- e(t))Cv* {-p) = -ECv*{-p), (25) 

which becomes the positive energy equation after factoring 
out the —1 and we have the charge conjugation relationship as 
before 



uip)^Cv*{-p). 



(26) 



Thus, this D2 alternative restores the invariance of the 
Dirac equation under charge conjugation transformations even 
in the presence of a vector potential. 



B. Ehrenfest Equations in D2 

As discussed in the introduction, the Ehrenfest equations 
in Dl for the classical Lorenz force depended on the sign of 
the energy of the wave function evaluated in the expectation 
values. A similar situation transpires for a spin in a constant 
magnetic field. 

The first of these is the calculation of the change in the mo- 
mentum in the presence of a vector potential and a scalar po- 
tential. In the Dirac Dl alternative the derived equation con- 
tains the sign of the en ergy of the eigenstates used to calculate 
the expectation values lUSll . The derivation is outlined in Ap- 
pendix A 



sgniE) — Umv)) 
at 



-e(E + v X B) 



(27) 



When the minimal coupling for alternative D2 is included, a 
second factor of sgn{E) multiplies the electric and magnetic 
fields on the right hand side and cancels out the factor on the 
left hand side. This recovers the classical Ehrenfest formula 
for the momentum 



-{{mv))^^e{E 



V X B). 



(28) 



In Appendix B a similar calculation for the spin dynamics 
of an electron in a magnetic field and factor of sgn{E) for the 
eigenstates forming the expectation value is also found in the 
same way as for the Lorentz force law in alternative Dl. In 



alternative D2 a second factor of sgn{E) is found on the right 
hand side and the appropriate Ehrenfest equation is recovered. 
So, D2 corrects the inconsistencies of the Ehrenfest equations 
inDl. 



C. Clioice of one electron basis states in D2 

In both alternatives the matrix elements of the Hamilto- 
nian can be determined relative to any chosen basis set. The 
usual approach is to begin with positive energy free electron 
plane waves. An equally valid basis set could be derived 
from eigenstates of hydrogenic atomic states or the eigen- 
states of any potential. The hydrogenic set of states would 
be more appropriate for situations involving a nucleus at the 
origin as found in heavy ion scattering and other atomic prob- 
lems. In the mass gap this hydrogenic basis will describe the 
bound states quite well. For energies outside of the mass gap 
these eigenstates, though derived from a nucleus at the origin 
could equally well describe scattering waves about the origin 
as well as other potentials. At large energies and large dis- 
tance from the nucleus the hydrogenic states are essentially 
coulomb scattering states which become approximately plane 
wave-like at very large distances. These basis states reflect the 
dominant influence of the charged nucleus. If there are other 
charges, as in scattering from other atoms, it may be desirable 
to distinguish between the internal potential which gives rise 
to the hydrogenic states and the external potential that could 
describe scattering in this basis based on the central nucleus. 
In both of the alternatives, the scattering states at high ener- 
gies and large distances will behave very much like the free 
electron bases. However, in D2 the atomic hydrogenic states 
will exhibit bound states in both the top and bottom halves of 
the mass gap because the energy specrum must be symmetric 
about zero. Electron and positron scattering results should be 
able to be described by any similar bases. The distinguishing 
characteristic will be the bound states in the gap. 



D. Bound States in tlie Mass Gap in Dl and D2 

When bound states in the mass gap are considered, there 
are very real differences between the two alternatives. In al- 
ternative Dl , there are only bound states in the positive half 
of the mass gap and no states in the lower half of the mass 
gap. In alternative D2 the positive energy bound states in 
the upper part of the mass gap are mapped by the charge con- 
jugation transformation into the lower half of the gap. The 
existence of these bound states in the lower half of the mass 
gap gives rise to a new set of excitations that do not exist in 
the Dl alternative. These negative energy bound states in the 
lower part of the gap must also be filled in the Dirac vacuum. 
Excitations out of these states will leave behind bound holes 
in the vacuum that are localized near the nucleus. The exis- 
tence of these states should give rise to excitations that were 
not expected in alternative Dl. It should be possible to distin- 
guish between the two alternative by examining the properties 
of these hole states. That is the goal of a later section in 



this paper which examines data from heavy ion scattering ex- 
periments mentioned above. D2 displays a different energy 
spectrum in the mass gap from Dl. 



E. Piece-wise Constant Potentials in Alternative Dl and D2 

The examination of constant potentials is important because 
any potential can be approximated mathematically by piece- 
wise constant potentials. In the following discussion the cal- 
culation will be carried out in one dimension for spinless par- 
ticles since it simplifies the mathematics. This reduces the 
matrices to 2 x 2's and spinors to 1 x 2's. 

Consider first a positive constant potential Vb in a certain 
region. In this region the wave function 



ip{z) = ae 



ikz 



Ul 
"2 



(29) 



will carry a current in the positive z direction. Now consider 
a comparison of this wave function and its energy spectrum E 
in both the Dl and D2 alternatives. 

In alternative Dl, the Dirac equation yields the following 
secular matrix equation 



Vq- E + ni k 

k Vq — E — m 



Ul 
U2 



0. 



It is easy to show that the energy eigenvalues are 



Ek =Va± 


Vm^ + fc2, 


and the wave functions are 


^P+{z)=ae'''' 


1 

k 




L m+Vm^+k^ -1 


and 




-fe -1 


Vfc-(^) = ae"=^ 


(rn + Vm^+k^) 
1 



(30) 



(31) 



(32) 



(33) 



The defining signs of the energies are relative to Vq and the 
actual sign of the energy will depend on the size of Vb . 

The "positive" energy states are found for E > Vo + m and 
the "negative" energy states are found for E < Vq — m and 
there are no states found in the mass gap VQ — m<E<V{) + 
m which is centered on Vq. As Vq becomes more positive, the 
"negative" (relative to Vq) states are pulled up in energy and 
can become degenerate with positive energy states in a nearby 
region, leading to Klein paradox effects. These effects were 
shown in Fig. 1 for the energy range {m, Vq — to) if Vb is 
large enough. Specifically, for a potential step of height Vb 
and extending infinitely far in one direction, the transmission 
coefficient T can be calculated analytically. If the potential 
step height Vq is less than the width of the mass gap, Vq < 2m, 
(in this section h — 1 and c = 1), the transmission coefficient 
T is zero for incident particle energies E in the range m < 
E < Vq + ra and becomes non-zero for E > Vq + ra. T 



1.2 pr 

1.0 - 

H 

I 0.8^ 

o 

^ 0.6- 

o 

I 0.4 - 
c 

H 

0.2 - 

0.0 t^ 



Vn-m 




Vo+m 




Energy E/m 

FIG. 1 . The one dimensional transmission Coefficient T (calculated 
in alternative Dl ) for a potential step of height Vo/m — 3 for < 
a; < oo plotted as a function of the incident particle energy E. The 
small peak between m < E <Vo — mis the Klein tunneling. 



approaches 1 as the incident energy becomes large. These 
results are very classical. 

If the potential step is larger so that Vb > 2m, the trans- 
mission coefficient develops a non-zero value between m and 
Vq — m. It is as if the barrier has become partially transpar- 
ent to particles in the range m < E < Vq — m. Mathemati- 
cally, this effect is due to the negative energy states (relative to 
Vq) being pulled up into degeneracy with the positive energy 
states of the incident particles. An example of the transmis- 
sion coefficient for this case is shown in Fig. 1, for the value 
of Vb/ni = 3. This situation of Klein tunneling presents a se- 
rious challenge for any theory of low energy particles because 
they would appear to tunnel into high and infinitely wide bar- 
riers with very little energy. This would seem to contradict 
widely observed classical phenomena. 

The situation for the D2 alternative is quite different. The 
secular matrix equation becomes 



sgn{E)VQ - E + m k 

k sgn{E)VQ — E — m 



and the eigenvalues satisfy 

sgn{E){\E\ - Vq) = ±Vm^ + k^. 
So, for the positive sign the eigenvalue is 

E+ = Vq + Vm^ + k^, 
while for the negative sign 



Ul 
U2 



-5,7 = -Vb - V^ 



fc2 



= 

(34) 

(35) 
(36) 
(37) 



and the coiTesponding wave functions are unchanged. 

In this D2 alternative the positive energy states are found 
for E > Vq + m and the negative energy states are found for 
E < —{Vq + m). There are no oscillatory states allowed in 
the intervening interval — (Vb + m) < E < Vq + m. Since 



there are no states in this interval, there is no possible cur- 
rent in this energy range. In contrast to the Dl alternative, 
in the D2 alternative increasing the strength of the potential 
simply widens the interval in which no states are allowed. It 
is not possible for "negative" energy states to ever be degen- 
erate with positive energy states in a neighboring region. The 
D2 alternative suppresses the Klein tunneling allowed in Dl. 



F. Klein Tunneling in D 1 and D2 

Consider the classic step potential of the Klein paradox 
(\/ = for z < and y = Vq for z > )in the D2 al- 
ternative. There are no states for z > in the energy interval 
m < E < Vq + m, and the transmission coefficient T must 
be zero in this interval. For energies E > Vq + m the trans- 
mission coefficient T is given by 



T 



Ar 



(1 



where 



'£;-(Vb + m) E + m 
E - {Vo - m)y E - m 



(38) 



(39) 



and E — y/m? + k^ is the incident energy. 

In the D2 alternative there is no Klein paradox. Increasing 
Vq simply expands the interval {m, Vq + m) where T = 0. 
If Fig. 1 were constructed in the D2 alternative, there would 
be no small peak between m and Vq — m and T = in that 
interval. There is no question about the potential being too 
strong. There is simply no current in the step in the energy 
interval {m, Vq + m). 

In this classic step problem it is possible to find evanescent 
waves whose energies are less than Vo + m. These evanescent 
waves do not carry currents, so their presence does not change 
the discussion about T vanishing in the interval {m, Vo+m). 

In the D2 alternative the evanescent wave function for z > 
0, for example, will have the following form 



V'(z) = ae '^^ X 
and the eigenvalues are found to be 



Wl 
W2 



E+ = Vo + Vm^ - 't^ 



and 



e: 



'Vq — \/ m? — K? 



(40) 



(41) 



(42) 



where < k < tti and for the positive energy state w\ ~ \ 
and Wl — in lira + \/m^ — k^). No states are allowed in 
the interval (— Vq, Vq). These evanescent states would enable 
tunneling in the energy interval (Vo, Vq + m), but not in the 
interval (— Vq, Vq)- 

It should be noted again that the energy spectrum of these 
states are symmetric about the zero energy. D2 does not allow 
the Klein tunneling effects which are in Dl. 



G. Attractive Potentials in Dl and D2 

In Dl a negative square well can have any depth. It is pos- 
sible to have bound states in the gap of such a potential at any 
depth in the mass gap. It was the question of what happens 
when the potential is deeper than the mass gap and bound 
states of this potential were to merge with the continuum of 
negative states below the mass gap that motivated the unsuc- 
cessful search for unstable vacua by Greiner and colleagues. 

In D2 the presence of sgn{E) gives rise to a surprising dif- 
ference from Dl. Consider a negative potential V{x) (for the 
moment restrict the potential to be one dimensional). Con- 
sider first a positive energy eigenstate of this system. This 
means that the sgn{E) factor multiplying the potential is pos- 
itive. For a negative potential the eigenvalues that are lower 
than m will be found in the mass gap (— m, m). How deep can 
the potential be and still give rise to a positive energy eigen- 
value? Clearly if the depth of the potential is less than m the 
bound states will be above zero and will clearly have positive 
energies. However, if a potential extends from +to down to 
some negative energy below zero, it is possible that the poten- 
tial could have a bound state that is below zero, that is nega- 
tive, but we have been considering only positive energies and 
this would lead to a contradiction. Thus, this reductio ad ab- 
surdum argument would imply that in alternative D2, a neg- 
ative potential cannot be deeper than —to or more precisely, 
cannot have bound states that are below zero energy. Because 
of the charge conjugation symmetry in D2, this would mean 
that the energy spectrum remains symmetric about zero and 
the eigenstates and eigenvalues are mapped from positive en- 
ergy onto negative energies. There cannot be any crossing 
through zero. 

From this perspective it is notable that the positive energy 
hydrogenic bound states of the Dirac equation are clearly in 
the interval (0, to,) and that the lowest possible state for a nu- 
clear charge Z is zero as Z ^ 1/ao, where a^ is the fine 
structure constant. 

In this D2 alternative a plot of the potential energies for 
both positive and negative energy states will have a reflection 
symmetry about zero energy and never cross over the zero 
energy. No attractive potential for negative energy states can 
extend above the zero energy level. There does not seem 
to be any such restriction on repulsive potentials since any 
eigenvalues would be larger than m. There does not seem to 
be any obvious restriction on the height of barriers. 

D2 seems to imply that no positive energy state can cross 
the zero energy line as is allowed in Dl. This property ex- 
plains why the search for an unstable vacuum was not suc- 
cessful. The ultimate question is whether nature actually ex- 
hibits distinctive behavior that confirms the D2 alternative as 
the correct alternative. 



H. Hole Theory in Dl and D2 

Let us now consider Dirac's argument for positrons in the 
hole theory in both Dl and D2. First, it must be re-iterated 
that in the free electron case there is explicitly no charge in 



the Hamiltonian. Charge only enters the description in the 
presence of a vector potential. Charge enters into the hole 
theory in an independent fashion through an argument for the 
change in charge of the vacuum in the presence of a hole. 

In the absence of any potential both Dl and D2 are the 
same. The Dirac vacuum with all of the negative energy states 
filled is necessary to avoid transitions to unbounded negative 
energies. In both alternatives a hole in the free electron vac- 
uum at an energy of — -Ep, momentum —p , and spin — s is 
clearly interpreted by the change in the energy of the many 
electron state as a positive energy state Ep_ with momentum p 
and spin s. The charge associated with a hole is imputed by 
observing the change in the charge of the vacuum plus hole 
in comparison with the vacuum. The charge in not carried 
by the wave function. The charge is the opposite of whatever 
the charge was of the original electrons whose wave functions 
were solutions of the free electron Dirac equation. 

D2 describes positron wave functions using charge conju- 
gation in the same way as does Dl. However, the positron 
wave function in D2 is identical with some positive energy 
electron wave function in contrast to Dl. The charge is not 
carried in the wave function, but is ascribed by the change in 
the vacuum as in _D1. 

In D 1 there is a long history of ascribing of the charge of the 
vector potential in the Dirac equation under the charge conju- 
gation transformation as evidence that the positron wave func- 
tion represents an opposite charge from the electron. From the 
perspective of D2 this sign change is simply the breakdown 
of charge conjugation invariance in Dl. In D2 the positron 
wave function must necessarily be identical to a positive en- 
ergy electron state just as is true for the free electrons. The 
positron charge is not carried inherently in the wave function 
and must be imputed separately. 



I. Current Density and Zitterbewegung In Dl 

A central argument for the Dl alternative has been the 
discussioni 19] of the probability current density for free elec- 
trons. The standard derivation for the continuity equation of 
the probability density yields for the probability current den- 
sity 



Ji = ca,. 



(43) 



where c is the velocity of light and a^ is one of the Dirac 
matrices. 

This result has long been regarded as peculiar since in the 
classical limit we expect the current density to be carried by 
the momentum pi/ra. A long standing puzzle related to this 
peculiarity is the fact that in the absence of a potential the mo- 
mentum is a constant of the motion. However, this probability 
current density does not commute with the Hamiltonian and is 
not a constant of the motion in the absence of interaction po- 
tentials. This is the third Ehrenfest equation in I?l which does 
not obey the classical limit. 

In the Dl alternative using a wave packet constructed out 



of positive energy states 

(44) 
it is straight forward to evaluate the expectation value of the 
probability current density and show that 



{Ji)+ = {ca,)+ = 



d^p 



2 2 

^—^2^\bip,sr^C-^)^ 



(45) 

where the subscript + on the averages indicates only the pos- 
itive energy states are used. 

At this point a crucial argument has been made in the Dl 
alternative. It is noted that the eigenvalues of the cai matrices 
are ±c, corresponding to positive and negative energies, and 
if the expectation value of the probability current density is to 
be calculated using the eigenvectors of the a; matrices it will 
require the inclusion of negative energy states in order to ob- 
tain a velocity associated with the probability current density 
less than the velocity of light. 

When the current density is evaluated using such a linear 
combination of states. 



V'(x) 



(2 



+d*{p,s)v{p,s)e'P-^X46) 

the resultfl91 is a weighting of the momentum by the ampli- 
tudes |6(p, s)p and \d{p,s)\'^, and also the famous zitterbe- 
wegung terms that exhibit frequencies on the order of the rest 
mass ±2mc^. These high frequency terms have been puz- 
zling and have led to a number of intuitive arguments about 
confinement and the excitation of electron-positron pairs in 
the presence of static potentials. 

The amplitude of the zitterbewegung in the expectation 
value of the probability current density is proportional to the 
amplitude d*{p, s) of the negative energy states in the wave 
function. The next series of Dl arguments are guided by an 
example wave function. A typical wave function 1 19il at zero 
time is 



V'(a;) = {nd'^y 



-3/4, 



'wi{0) 



(47) 



where d is the "confinement width" of the wave function and 
^1(0) is the spin up spinor for an electron with zero momen- 
tum and positive rest energy. 

The amplitude of a positive energy plane wave in this state 
is 



6(p,s) = .4(^)3/V 






and the amplitude for a negative energy state is 
d^ 



d*{p,s)=A{^r/'e- 



where ^ is a normalization factor 



pc 



Er, 



mc^ 



(48) 



(49) 
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The relative fraction of the negative energy states in this 
wave function is small unless pc « mc^ or greater At the 
same time the confinement of the wave function to a region of 
size d implies the momentum must be p < h/d, or that the 
confinement is comparable to the Compton wavelength 



d: 



n 

mc 



(50) 



From these arguments in Dl there have evolved two con- 
clusions. If the confinement of a wave function is on the 
order of the Compton wavelength or smaller, negative en- 
ergy states will be significantly probable. Because of the 
connection between negative energy states and positrons this 
would mean that static, confining potentials should generate 
electron-positron pairs under the conditions of close confine- 
ment. 

The origin of these arguments in Dl is the fact that the 
probability current density is not proportional to the momen- 
tum operator and to get a speed associated with the probability 
current density to be less than the velocity of light, negative 
energy states must be mixed into any wave function. The 
relevant question, besides the question of whether all of the 
inconsistencies of I?l are acceptable, is whether the probabil- 
ity current density is what transports the charge density. This 
question is resolved in D2 in the next section. 

These ideas are deeply engrained in those of us who have 
learned quantum mechanics in alternative Dl. In particular, 
we have been forced into these outcomes by the form of the 
probability current density. Let us now examine the same 
picture from the point of view of alternative D2. 



J. Charge Current Density in D2 

First, it must be noted that both alternatives Dl and D2 
will have the same probability current density for free elec- 
trons. For free electrons there is no distinction between the 
two alternatives. However, there is also explicitly no charge 
in the free electron Hamiltonian. To examine the charge cur- 
rent density it is necessary to examine the Hamiltonian that 
couples the electrons to a vector potential. 

The D2 alternative differs primarily in the treatment of the 
electron equations in the presence of a vector potential yl^ = 
{^,Ai). For a Hamiltonian that couples vector potentials to 
the free electrons, the charge current density jf^ is most easily 
derived by 



j: 



Q 



dH 



(51) 



Since the Hamiltonian in alternative D2 is 



H ~ ca ■ p + jSmc^ — e sgn{E)cx ■ A + e sgn{E)^, (52) 
the charge current density is 
1 



where in the second equality the charge current density op- 
erator has been symmetrized as needed for the free electron 
Hamiltonian. For free electrons with momentum p , we have 
that sgn{E) can be represented by a projection operator 



sgn{E) 



ca ■ p + jSmc^ 

Er, 



(54) 



Substituting the expression for the free electron projection op- 
erators into the expression for J^ we obtain 

J^ = e{l/2){ca ■ pai + caiOc ■ p + mc{/3, 0:^}) (55) 

and using the anti-commutation relations for the ai and (3 ma- 
trices 



{ai.aj} = 25i, 



{a„/3} = 



(56) 



(57) 



we obtain the current density operator proportional to a com- 
ponent of the momentum operator 



J^ 



hec 



Ep' 



(58) 



J> 



esgn{E)ai — -e {sgn{E)ai + aiSgn{E)), (53) 



the matrix li is the 4x4 identity matrix and pi is the i-th com- 
ponent of the momentum operator 

Thus, in alternative D2 the charge current density is in ex- 
cellent accord with classical expectations. It is proportional to 
the momentum, is a constant of the motion in the absence of a 
potential, thus satisfying the Ehrenfest equations which were 
violated in alternative Dl. Because the charge current density 
is proportional to the momentum there is no zitterbewegung ! 
There is no paradox as in Dl where the probability current 
density operator has eigenvalues which force the inclusion of 
negative energy states in wave packets. Expectations of the 
current density can be determined completely by positive en- 
ergy wave packets as above. There is no necessity for the 
inclusion of negative energy states in wave function packets 
as was argued in Dl. . In alternative D2 the negative energy 
states are derived from the positive energy states by the charge 
conjugation transformation and only the positive energy states 
are necessary for completeness. 

The arguments that confinement by static potentials neces- 
sarily induces particle anti-particle states in so far as it was 
based on arguments of eigenvalues of ca no longer has any 
imperative in alternative D2. In D2 the charge current den- 
sity, which is what describes the movement of the charge, is 
proportional to the momentum and there is no zitterbewegung 
or generation of positrons by confinement in static potentials 
as argued in Dl. 



K. Completeness In D2 

It is often argued in Dl that positive energy states are not 
complete. This argument relies completely on the apparent 
need in Dl for negative energy eigenvalues of the a matrices. 



10 



As has been discussed above there is no need for this argument 
when considering the charge current density. In D2 the stan- 
dard mathematical proofs based on the Fourier Theorem n24l 
can be used to prove that the positive energy states are com- 
plete in the Hilbert space. Since the negative energy states 
are connected to the positive energy states by the charge con- 
jugation transformation, they are not linearly independent of 
the positive energy states. This means that the positive energy 
states are complete and are the only set of states needed. The 
inferences in Dl based on the probability current density and 
the eigenfunctions of the a matrix are irrelevant to the charge 
current density and thus the movement of charge in D2. The 
arguments on confining potentials generating e^e+ pairs does 
not arise in D2. In D2 the positive energy states are complete 
in contrast to the conclusions in Dl. 



Transformations and Feynman/Stuckleberg Theory in Dl 
andD2 



In both of the alternatives Dl and D2 there are three 
transformations that can be formulated in the free electron 
basis. These transformations are well defined in several 
references lllSl [l9l l25ll . They are: C charge conjugation, 
T time reversal, P parity transformations. If we write the 
Hamiltonian in two parts: 



Hq = a ■ p + f3m 



(59) 



and 



Hi = -ea. • A + e$, (60) 

it is easy to show the following transformations. 

CHo{p)C-^ = -Ho{-p), (61) 



PHo{p)P-' = Hoi-p), 



THoip)T-' = -Hoi-p). 



(62) 



(63) 



The first of these shows that the free electron Hamiltonian has 
a reflection symmetry about zero energy and is charge conju- 
gation invariant. By exactly the same process for free elec- 
tron states, we can show 



CsgnE{p)C- 



-sgnE{-p), 



PsgnE(jp)P ^ = sgnE{—p) 



TsgnE{p)T ^ — —sgnE{—p) 



(64) 



(65) 



(66) 



If we apply these transformations to the interaction part of 
the Hamiltonian, we obtain 



CHiC- 



Hi, 



PHiP-^ = Hi, 



tHit-^ =Hi. 



(67) 



(68) 



(69) 



In alternative D2, we have charge conjugation invariance 
and have the following transformations. 

CiHo{p)+sgnE{p)Hi)C-' = -{Ho{-p)+sgnE{-p)Hi, 

(70) 

T{Hoip)+sgnE{p)Hi)T-^ = -{Ho{-p)+sgnE{-p)Hi). 

(71) 
Negative energy wave functions are transformed into positive 
energy wave functions and vice versa. The charge conjuga- 
tion wave function ^cOf a negative energy state ?/;(_) is given 
by 



V-c = cv-, 



(-) 



(72) 



is identical with a positive energy state with the appropriate 
momentum and spin and define the positron wavefunction. 
Just as in alternative Dl, because the time reversal transfor- 
mation gives rise to the same transformation properties as for 
the charge conjugation, the time reversed wave function will 
necessarily be the same (within multiplication by a constant 
matrix) as the charge conjugation transformed wave function 
except for a phase factor The Wigner time reversed wave 
function "ii pcT 1.19.1 is an explicit representation of this con- 
nection between positron wave functions and a matrix multi- 
plied by the wave function of a negative energy electron mov- 
ing backward in space time and is valid in D2 as well as Dl. 
In alternative Dl charge conjugation invariance is broken 
and we have 



C{Ho{p) + Hi)C-^ = -(Hoi-p) - Hi). 



(73) 



In much of the literature about alternative Dl this breakdown 
of the charge conjugation invariance has been interpreted as 
evidence for the sign change of the positron. This has per- 
sisted even though the positron charge was determined by the 
vacuum neutrality condition. From the perspective of the D2 
alternative, the change of sign is just the breakdown of charge 
conjugation invariance. 

In D2 it is possible to demonstrate the apparent charge dif- 
ference by examining the negative energy wave function and 
its time reversal transform. For a positive energy, the wave 
functions must obey 

'^~m = (" • (P " ^^^) +Prn + e$)V'(+) ■ (74) 

For a negative energy the wave function must obey an equa- 
tion that looks like an oppositely charge particle, except for 
having a negative energy. 
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Dl Alternative 



D2 Alternative 



ih 



dt 



{a- {p + eA) + l3m - e$)V'(_) 



(75) 



The time inversion transformation identifies the transforma- 
tion of this negative energy wave function as a positive energy 
wave function traveling backward in time. 



M. Feynman Diagrams and Perturbations in Dl and D2 

Because in D2 the negative energy states are determined 
by mapping from the positive energy states using the charge 
conjugation transformation and the positron amplitude has 
been identified with negative energy wave functions moving 
backward in space time, the application of the whole dia- 
grammatic structure of the Feynman Stuckleberg perturbation 
structurel 19J should be exactly the same in Dl and D2 when 
applied on positive energy states. The negative energy results 
could then be obtained by the charge conjugation transforma- 
tion. 



IV. HEAVY ION SCATTERING AS EVIDENCE FOR 
ALTERNATIVE D2 

In this section the experimental energies of the sharply de- 
fined e~e~^ pairs seen in heavy ion scattering will be com- 
pared with predictions of the D2 alternative. The properties 
of D2 delineated in the previous sections already predict that 
no evidence of an unstable vacuum would be found because 
positive and negative energy eigenvalues are separated by the 
zero energy. 



A. Historical Introduction 

P. A. M. Dirac's tlS] proposal for a "vacuum" in which 
all of the negative energy electron states were filled has led 
to a detailed understanding of free positrons, including their 
production and annihilation! 19]. However, the extension of 
his insights to bound, hydrogenic-like atomic states of the 
Coulomb potential has resulted in a strange story of contra- 
dictory experimental and theoretical results. Using a widely 
held version of alternative Dl and vacuum for bound states 
P] , Greiner and numerous colleagues [ 1 1 , 1 2 1 proposed a series 
of compelling effects to be seen in the high electric fields of 
heavy ion scattering experiments. 

These predictions and experiments were the motivation for 
three international proceedings: Landstein[4| MarateaH and 
Cai-gese|6]. 

The 1990 summary lecture [261 by Miiller, for the last in- 
ternational meeting [6], assembled a table of e^e+ pair sum 
energies and noted that "the phenomenology of data was so 
complex that they do not fit into any simple scheme.." and 
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FIG. 2. Representation of the types of electron positron excitations 
available in Alternatives Dl and D2. The two excitations on the 
left are available in both Dl and D2, but the excitation on the right 
represents a bound e^e^ pair that is only possible in alternative D2. 



little hope was offered of finding such a scheme. An intro- 
ductory essay for this same conference by J. Rafelski lBTII ru- 
minated that striking experimental observations were needed 
showing the breakdown of the vacuum in order to be sure that 
their understanding not be "Ptolemean", that is, making a fun- 
damental mistake in a basic notion, leading to more and more 
complex descriptions. None of the theories in this meetinglg] 
were able to make sense of the carefully crafted experimental 
data. 



B. Hydrogenic bound e e states in the mass gap 

The charge conjugation invariance in D2 gives rise to a very 
different structure of states in the mass gap. The positive en- 
ergy states for the hydrogenic atom are the same in both Dl 
and D2. In D2 the charge conjugation invariance maps the 
positive energy bound states to the lower half of the mass gap. 
Shown in Fig. 2 is a representation of the mass energy gap 
of the Dirac equation in both the Dl and D2 alternatives. In 
each alternative the Dirac vacuum will be formed by assum- 
ing all of the negative energy states are occupied by electrons. 
In the D2 alternative that would include the negative energy 
bound states in the lower part of the gap. Of course, in the Dl 
alternative there are no bound states in the bottom of the gap. 

The excitations on the left side of Figure 2 are present in 
both Dl and D2 because they involve holes below the mass 
gap. The left most excitation creates a free e^e+ pair. The 
second excitation creates an excited atomic state and a free 
positron. 

The right side of Fig. 2 is only found in the D2 alternative 
and illustrates an excitation from an occupied negative energy 
bound state to an (unoccupied) positive energy bound state. 
This state will be described as a bound e~e+ pair state. These 
bound e~e^ pairs within an ion constitute a new excited state 
of the system. The movement of these bound e^e+ pairs 
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is carried by the ions. Similarly, vertical excitations during 
ion-ion collisions from the filled vacuum states to an empty 
(ionized) positive energy state would give rise to a metastable, 
bound e^e+ pair state on one or both of the scattered ions. 

In the following we will argue that these excitations are re- 
sponsible for the narrowly defined free e^e+ pair energies 
observed in heavy ion scattering. For that argument it is first 
necessary to determine the excitation energies for these bound 
e^e+ pairs. 

In the D2 hydrogenic approximation the energy levels of 
the e^e+ bound pairs will be the difference between the posi- 
tive energy bound hydrogenic state E+ (n, j) and the negative 
energy bound state E^ {n,j) where 

E±{nJ) = 



±m 



l + ( 



Zuq 



n - U + i) + ^{j + 1)2 - {Zaof 



where Zuq is the product of the charge on the nucleus and 
the fine structure constant. In the hydrogenic approximation 
the bound e^e"*" pair energy difference would be 

Ae{n,j;n',f)^Ae{S;S') 

^E+(n,j)-E^{n',f) 
^E+{n,j)+E+{n',f) (77) 

where in this study the quantum numbers of the negative en- 
ergy bound states are primed and the atomic states will be 
labeled by atomic shell structure notation S = K,L1,L2, .... 
To get a better idea of the consequences of these new bound 
electron positron pairs in the D2 alternative, consider a Pb 
atom as an illustrative example. This means that Z = 82 and 
we can calculate the lowest three levels: 



E±{lSi/2) = E±{K) = ±409.4 KeV, 



E±{2S 



l/2j 



E±{Li) ^ ±484.9 KeV, 



E±{2R 



3/2; 



E±{L2 



±487.6 KeV. 



(78) 



(79) 



(80) 



In the comparisons with experiments below, we will use 
the K,Li,L2 notation for these positive energy states and 
K',L[,L2 for negative energy states. The labels for the 
metastable bound e^e+ pairs can be labeled by the atom, 
the positive energy and the negative energy, for example, 
Pb : K ^ K' . This designation would represent the decay of 
the metastable state consisting of a positive energy electron in 
state K and a negative energy hole from state K' . The lowest 
five different possible energies of different metastable states 
using the lowest three atomic states for Pb are shown in the 
table below. 



Transition Pb 


Energy Ae{S; S') 


Pb:K^K' 


818.8 Xey 


Pb:K^ LV 


894.3 KeV 


Pb: K ^ L2' 


897.0 XeV 


Pb:Ll^ LV 


969.8 KeV 


Pb:L2^ L2' 


97^.2 KeV 



If these metastable states can be created in heavy ion colli- 
sions, then there should be a variety of discrete energies avail- 
able for both beam and target atoms as the metastable states 
decay. 

In the D2 alternative when the two nuclei are close together 
during a collision, the positive and the negative energy states 
will be pulled toward zero. Since there is so much energy 
around during a collision, it is expected that many low lying 
positive energy states will be empty and electrons can be ex- 
cited out of the vacuum into these positive energy states leav- 
ing behind a hole in the negative energy states. When the nu- 
clei are near to their closest approach the difference between 
the positive and negative energy levels will be its smallest and 
the probability of excitation correspondingly greater As the 
nuclei separate in the scattering event one or both of the ions 
could be carrying the bound metastable e^e+ states that are 



(76) present in the D2 alternative. 



C. Decay Modes for bound e e+ metastable states 

Experimentally, one could probe the various decay modes 
of these metastable excited states and identify the various 
states by measuring the energy of the decay modes. In order 
to do this we first need to examine some of the decay modes. 

These metastable bound e^e+ states can decay in a num- 
ber of ways. In order to describe these different channels, 
we adopt a modified Feynman diagram description in which 
we denote the ion lines as a broader line if the ion contains 
a bound e~e+ pair and a slightly narrower line if there is no 
such excitation in the ion. The electron and positron lines will 
be thin lines. 

The strongest channel for decay of these metastable states 
would be the annihilation of the bound e^e+ pair state and 
the emission of a single photon. Because the ion can take up 
the recoil momentum, a single photon decay is possible. This 
is displayed in part A of Fig. 3. In these diagrams the time 
moves from left to right and the change in the width of the 
ion line represents the initial presence and then absence of the 
metastable state. This mode of metastable state decay can be 
either spontaneous, stimulated, or induced by other interac- 
tions. As will be discussed later, at high beam currents, in 
the presence of many ions, e^e+ pairs, and quickly changing 
fields, it might be expected that the metastable state lifetime 
is shortened. 

An example of other decay mechanisms can include the cre- 
ation of free e^e+ pairs. A lowest order process is shown in 
part B of Fig. 3. One would expect this process to give a 
broad range of pair energies starting at the 2m threshold. This 
process would not be expected to be a source of the sharply 
defined free e^e+ pair energy distribution. This process will 
require quite a large energy change in the ion kinetic energy 
to reach the pair production threshold. 

Another way in which the scattered, metastable ion could 
create a free e^e+ pair would be to join with an available 
photon. The diagram for this process is found in part C of 
Fig. 4. Here again it is hard to see how this process could 
give rise to a narrow kinetic energy pair peak. This process is 
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FIG. 3. Relaxation mechanisms for ion with bound e~ e^ pair (thick 
line) and without (thinner line). A: emission of photon upon anni- 
hilation of bound e~ e^ pair. B: emission of photon and creation of 
free e~ e^ upon annihilation of bound e~ e^ pair. 





line) and without (thinner line). C: emission of a free e" e^ pair in 
conjunction with absorption of a photon and annihilation of bound 
e~ e^ pair. D: emission of a free e^e^ pair in conjunction with two 
virtual photons and annihilation of a bound e^e'^ pair. 



a twisted form of the Bremsstrahlung scattering in a coulomb 
field, but with the ion recoil and bound e^e+ decay. Figure 
4D shows a process in which the pair energy must come from 
the change in the ion kinetic energy and the bound pair energy. 
It is not expected produce narrow free pair lines. 

There is at least one process that combines a possibly 
sharply delineated total energy with both the threshold 2m and 
the metastable bound state energies Ae(S', S'). It is of sixth 
order and one of the 16 diagrams is shown in Fig. 5. There 
are two vertices at which the excitation energy Ae(S', S") can 
be transferred and for each of these there are 8 different way s 
the virtual photons can be arranged. 



D. Evaluation of the cross section as a function of energy 

The evaluation of the amplitude and the cross section for 
this diagram raises a number of specific issues. The center 
of mass of the scattered ion can be treated as a Dirac particle. 
The electron/positron states in the diagram in lowest order can 
be treated as free particles. The internal degrees of freedom of 
the atom must be treated in terms of hydrogenic solutions of 
the Dirac equation and in calculating the cross section energy 
projection operators for the hydrogenic Dirac equation must 
be included in the product of operators. The full evaluation 
of the diagrams will not be included in this paper Here we are 
only going to discuss the general structure of this diagram and 




FIG. 5. Relaxation mechanisms for ion with bound e~e pair (thick 
line) and without (thinner line). Annihilation of a free e~e^ pair 
emission of a free e^e+ pair in conjunction with three virtual pho- 
tons and annihilation of bound e^e^ pair. 



examine the energy at which it might be possible to observe 
sharp lines in electron positron energies. 

The first such consideration is that we expect the annihila- 
tion of the initial free electron positron pair to be most impor- 
tant for very low energies. It is important to remember that 
in the limit of small relative velocities Vrei of the electron and 
positron pair the cross section for pair annihilation blows up 
at very small Vrei and has the form lll9ll 



Const 

Vrei 



(81) 



We should expect that the cross section for this process 
should be able to be written as 



fe+d 



E 



wr 



'|2 



Vrei 



-5{Ef ~ E.,), 



(82) 



where |Af' p is the square of the invariant amplitude and all 
of the other factors involved in evaluation of this amplitude, 
Ef is the final state energy, and Ei is the initial state of the 
process. The summation indicates the sums over interior co- 
ordinates. 

The final state energy includes the kinetic energy of the final 
free electron pair r+ + T_ , the final kinetic energy of the ion 
KEf, and the final internal energy of the atom E{Sf) {Sf is 
the designation of the atomic electronic states), and m is the 
electron rest mass. 



E 



f 



T++T^ + 2mc^ + KEf + E{Sf). 



(83) 



The initial energy will depend on the initial kinetic energy of 
the ion KEi , the initial internal energy of the atom E{Si), 
and the kinetic energy of the initial free electron-positron pair. 



14 



This kinetic energy can be written for the pair as a center of 
mass contribution Kcm and the relative kinetic energy Krei 



E, = Kern + Krci + 2mc^ + K E, + E{Si 



(84) 



Since we are looking at low energies we consider an average 
over the relative velocities of the free e^e+ states. In the low 
energy range of the initial annihilating electron positron pair, 
the relative kinetic energy is classical, Krei — \M{vrei)'^, so 
that an average over the initial velocities can be written as a 
sum over the initial relative kinetic energies Krei ■ 

a - / D{Krel)^J==S{Ef - E,)dKrel, (85) 

where D{Krei) is proportional to the density of states of the 
low energy electron positron pairs present during and shortly 
after the closest approach of the ionic collision. Since K^ei 
must be positive and 

Ef-E, = r++T_-Ae(5,; Sf)+^KE-K„a-Krei (86) 

where 

^e{Sr-Sf)^E{S,)-E{Sf), (87) 

and 



and 



we may write 



l^KE ^ KEf - KE, 



Se ^ AKE - Kc 



D{x)\M'\'^Q{x) 



where 



= T++T^-Ae{S^■Sf) + 5e 



(88) 



(89) 



(90) 



(91) 



and Q{x) is the unit step function. This cross section as a 
function of the final free electron positron kinetic energy T+ + 
r_ has an infinitely high and sharp peak at 



r+ + r_ = Ae(S'j; Sf) - AKE + K^ 



(92) 



To estimate the lowest possible energy to compare with ex- 
periments we will assume that Kcm is negligable. This as- 
sumption restricts our results to values close to the threshold. 
The pairs with non-negligable values will be distributed over 
a range and will thus not contribute to a sharp peak. 



E. Free e e+ Pairs in the Ion Center of Mass Frame 

The expression for the kinetic energy in the center of mass 
frame of the system will be 



In what follows we do neglect the kinetic energy Kcm of this 
initial pair because we seek the lowest energy situation. 

In order to compute the kinetic energy of the free e^e+ pair 
in the lab, first we must calculate the kinetic energy of the ion 
in the center of mass frame. Let vj be the magnitude of the 
lab frame velocity of the ion while it is in an excited state. In 
the laboratory the ion kinetic energy would be \Mv'j ,then the 
kinetic energy in the center of mass will be ^Mvj/^j where 
7/ = l/-\/l — /?|, M is the mass of the ion, /3/ ~ vi/c, and 
c is the velocity of light. 

Initially in the center of mass frame the energy of the ion is 



1 



E[ = 2mc^ + Mc^ + -Mvjhi + Ae(5; S") 



(94) 



and the momentum is zero. 

To find the lowest possible energy of the emitted free e~e^ 
pair, we consider the momentum change MAvj of the ion to 
be in the same direction as the initial velocity and the momen- 
tum of the free e^e+ pair will be in the opposite direction and 
the pair will have an opening angle of 29^ ■ The momen- 
tum equation in the CM frame for the ion before and after the 
emission of the free e^e+ pair is 



= MAvi - 27emc,3eCos(6'e). 



(95) 



The final energy after the emission of the free e e+ pair is 

E'^ ^ 2mc^ + -M{vi + Avif + (7^ - l)2mc^ . (96) 

Neglecting the term with (Aw/)^ which is very small, we have 
the equation 

Ae(5'; S') = (7e - l)2mc^ + ^^^-^cpi. (97) 

7/ 

Using the conservation of momentum equation to solve for 
MAvi and dividing by 2?7ic^ results in 



AejS; S') 
2mc^ 



where 



= 7e - 1 + lefieR, 



Pi 



R— — cos(6'e). 

II 

This results in an equation for 7e of the form 

Ae(S'; S' 



(98) 



(99) 



2?TiC^ 

This equation has two solutions for 7, 

/-, A£(S;S') x 
^^ 1 - i?2 



1 - 7e = R^Jil - 1. (100) 



TZ' = Ae(5; S') - AKE^^n + K, 



(93) 



±-^J^i(^(2 + M^) + ^..(:oi) 



1 - i?2 V 2mc2 



2mc^ 
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F. Transforming to the laboratory Frame 

In the center of mass frame for the ion, the e^e+ pair has 
the energy i<^' _ + and momentum P'e-e+ given by 



K-^+ = (7e - l)2mc2 + 2mc2 



and 



P'e-e+ = — 2mc^7e/3e COs( 



(102) 



(103) 



The pair energy T 
transformation 



(lab) 



in the lab will be given by the Lorentz 



2mc^ + T 



(lab) 



-n(E'^-^++PiPi,-,+ ). 



(104) 



From the center of mass energy equation we have 



(7e - l)2mc2 = Ae(5': S") - 2mc^-fe(3e—cos{9e) (105) 

7/ 

and we obtain for the total pair kinetic energy in the laboratory 
frame 



T 



(lab) 



= (7/-l)2mc2 + 7,(Ae(^;5') 



-2mc' 



1 + 7J 
7/ 



v/7|-l/S/cos(0e)). (106) 



This is the expression that should be compared with the exper- 
imental peak locations. It depends on the excited state energy 
of the bound e^e+ pairs, the initial velocity of the scattered 
ion, and the opening angle 26*2 of the free e^e+ pair in the 
center of mass frame of the ion. 

The excited bound pairs on the ion are created during the 
close approach of the scattering event and the emmision of the 
free e^e+ pairs must occur after the ion has been scattered. 
Because most of the spectrometers are constructed to measure 
the e~e"'" pairs perpendicular to the initial beam direction it 
suggests that only the ions that have scattered through small 
angles will emit e~e^ pairs that could be detected. Since 
there is no knowledge about the scattering angle of the ions 
in most of the experiments, we have here assumed that the 
velocity of the ions after scattering is approximately equal to 
the magnitude of the initial beam velocity for all of the ex- 
perimental situations. If the scattering angle of the ion was 
known it would be possible to solve for the scattered velocity 
in terms of the initial beam velocity and correct this assump- 
tion. Most of the experiments were conducted in a range close 
to X = 6 Mev per nucleon. Accordingly, in the comparison 
with experiment, we will use 



7/ = 1 + O.OOlx 



(107) 



and also determine /?/ from this value for x — 6. 

The two solutions for the je depend on the opening angle 
26'e of the free pair in the center of mass frame. This angle 
ranges from to tt/2. The value of 6*6 = is unphysical 
since the two particles would be moving in the same direction 
on top of each other The value 9^ = 7r/2 would have the 
electron and positron moving in opposite directions and that 
would allow no recoil of the ion in the center of mass frame. 
The physical values must be in between. 



G. Comparing Experimental Peak Energies and Free e e"* 
Lab Energies 



The experimental data has been organized in two tables. 
Table 1 contains data from spectrometers that measured both 
the electron and positron energies. Table 2 contains data from 
earlier experiments in which only the positron was observed. 
We have first calculated all possible total pair kinetic energies 
for each ion (beam and target ) of each experiment. We ini- 
tially chose the value of the transition at 6e = 7r/4 which was 
closest to the experimental peak energy. This identified the 
possible transitions that could have contributed to the peak. 
We then solved for the exact value of 9^ that would have re- 
produced that peak. In some cases there were more than one 
transition possible and we have listed them all. If a transition 
had a very different value of 9^ it is likely not to be the cor- 
rect transition. Since there were two different solutions of the 
equation for 7e each transition will be labeled with the sign 
of the solution and the atomic shells. For example, one of 
the smallest transitions would be labeled (+) U : K ^ K' . 
Most of the experiments were well described by shells of 
K, LI, L2 for positive energy states and K' ^ LI', L2' for neg- 
ative energy states. For some of the nigher energy peaks it 
was necessary to use (M : n — 3,j — 1/2) and for the 
very highest transitions {Z : n Ri 50, j — 1/2). Most of 
the experiments are well described by the lowest transitions: 
K ^ K',K ^ L1',K -^ L2',L1 -^ L1',L2 -^ L2' and 
the peaks were arranged in the corresponding order The fact, 
that in the U + Pb experiment we identified some small peaks 
from the experimental data that would have never been re- 
ported by a good experimentalist, and these fit nicely with ap- 
propriate transitions, is gratifying. 

Each line of the tables includes, besides the experiment 
identifier, the following data: ( 1 ) the experimental peak value, 
(2) the theory peak value for 9^ — 7r/4, (3) the sign (±) for 
the 7e solution, (4) the transition A : X ^^ Y, and (5) the 
value of 9e which gives the experimental peak value. 

Every experimental line was identified with at least one 
plausible transition. Each of the tables is discussed in some 
detail in some remaining paragraphs. 

The remaining issue in the heavy ion scattering story his- 
torically is the experimental observation of the APEX col- 
laboration that as the beam current was increased, the sharp 
e^e+ peak slowly "melted" into the background. The con- 
clusion that the peaks were thus spurious and unexplainable 
is based on the assumption that the peaks represented a stable 
entity. However, as seen above, the bound e^e+ pair states 
are metastable and the stimulated decay rate would certainly 
increase as the number of scattered ions and other scattering 
products increase with beam current. Appendix C contains ar- 
guments that imply the induced decay of the metastable bound 
e^e+ states will increase with beam current. As the induced 
decay rate increases, there will be fewer free e^e+ pairs to be 
observed. This metastable state effect could explain why the 
APEX experiment failed to see a clear e^e+ energy peak as 
the beam current was increased. 
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TABLE I. A comparison of the experimental 
free e~ e^ pairs emitted from heavy ions and 
tions for specified ions. 



sum-kinetic energy of 
the theoretical transi- 



System,Ref. 


Obs. 


Theory(6le=45°) 


± 


Transitions 


Oe 


U+Pb[28] 


576 


521.4 


+ 


U:K^K' 


56.0° 






571 


+ 


Pb:K^K' 


46.4° 




680* 


677.6 


- 


Pb:K^Lr 


46.0° 






679.8 


- 


Pb:K^L2' 


45.4° 






680.9 


+ 


U:L1^L1' 


45.3° 






687.8 


+ 


U:L2^L2' 


43.9° 




734* 


740.8 


- 


Pb:Ll^Ll' 


43.5° 






745.3 


- 


Pb:L2^L2' 


42.4° 




787 


784.5 


- 


Pb:Z^Z' 


69.2° 




934* 


784.5 


- 


U:Z^Z' 


69.3° 


U+U[28] 


553 


563.5 


- 


U:K^K' 


42.0° 




634 


645.2 


- 


U:K^Lr 


42.2° 




634 


649 


- 


U:K^L2' 


41.0° 


Th+Th[2E] 


595 


574.5 


- 


Th:K^K' 


51.1° 






574.5 


+ 


Th:K^Ll' 


47.9° 




608 


607.5 


+ 


Th:K^Ll' 


43.8° 






610.8 


+ 


Th:K^L2' 


44.8° 


U+Th[3a] 


760±2 


763.3 


- 


Th:M^M' 


44.6° 






762.2 


- 


U:M^M' 


44.8° 


U+Ta[28J[30] 


630±8 


645.1 


- 


U:K^L1' 


40.9° 






649 


- 


U:K^L2' 


39.8° 




746 


750.9 


- 


TaiLl^Ll' 


44.2° 




805±8 


784.5 


- 


Ta:Z^Z' 


50.0° 



TABLE II. A comparison of experimentally measured free positron 
energies and predicted energies from designated transitions for spe- 
cific ions. The theoretical positron energy was assumed to be half of 
the total kinetic energy of the pair. 



H. Discussion of Table I 



System,Ref. 


Obs.£+ 


E+(e^ = 45°) 


± 


Transition 


Oe 


U-l-Cm(epos) 


328±9 


322.6,324.5 


- 


U:K^L1',L2' 


48,47 




328±9 


335.5 


-1- 


Cm:Ll-j.Ll' 


42.6 


U-l-Cm(epos) 


445±12 


392.2 


- 


Cm:Z^Z' 


67 


Th-l-Cm(epos) 


354±10 


365.4 


- 


Th:Ll^Ll' 


39.3 




354±10 


358.7 


- 


Cm:Ll->Ll' 


42.9 


Th-l-Cm(epos) 


367±9 


363.9 


- 


Cm:L2^L2' 


46.9 


Th-l-Cm(epos) 


420±11 


392.2 


- 


Cm,Th:Z->Z' 


57.6 


U-l-U(orange) 


283 


281.7 


- 


U:K^K' 


46.1 


U-^U(epos) 


354 


363.6 


- 


U:K^K' 


40.2 


Th-l-U(epos) 


367 


368.5 


- 


Th, U:L2^L2' 


45.2 


U-l-Th(orange) 


291 


281.8 


- 


U:K^K' 


50.7 




291 


287.3 


- 


Th:K^K' 


47.6 


U-l-Th(epos) 


354 


344.9 


- 


Th:L2^L2' 


48.5 


U-l-Th(epos) 


459 


392.2 


- 


U, Th:Z^Z' 


71.9 


Th-l-Th(epos) 


314 


326.1,327.8 


- 


Th:K^Ll',L2' 


38,37 


Th-l-Ta(epos) 


367 


368.5 


- 


Th:L2^L2' 


44.6 




367 


375.5 


- 


Ta:Ll^Ll' 


40.9 


U-l-Au(orange) 


261 


260.7 


-1- 


U:K^K' 


45.5 




261 


281.8 


- 


U:K^K' 


30.4 


U-l-Au(orange) 


327 


320.3,321.2 


-1- 


Au:K^Ll',L2' 


48,47 




327 


322.6,324.5 


- 


U:K^L1',L2' 


48,47 


Pb-l-Pb(orange) 


331 


338.8,339.9 


- 


Pb:K^Ll',L2' 


41,40 


U-l-Ta(orange) 


302 


304.2 


-1- 


Ta:K^K' 


44.3 




302 


300.3,302.2 


-1- 


Ta:U^Ll',L2' 


46,45 



Table I contains the measured peak locations for a number 
of the experiments which have been published in the literature. 
A peak energy with a * represents a less than significant peak 
identified by the authors of this theoretical study and would 
probably not be regarded as significant by the original authors 
of the experiment. Each line of the table includes, besides the 
experiment identifier, the following data: (1) the experimen- 
tal peak value, (2) the theory peak value for Be = 7r/4, (3) the 
sign (±) for the 7e solution, (4) the transition A : X ^ Y, and 
(5) the value of 6'e which gives the experimental peak value. It 
is satisfying that the less significant peaks in the Pb-nPb exper- 
iments are fit so well to the low lying level of the Pb atoms. 



I. Discussion of Table 11 

The experimental efforts to study these effects were car- 
ried out principally by two groups whose spectrometers were 
given the labels: EPOS and ORANGE. In the early manifes- 
tations of these spectrometers only the positron energy could 
be determined. Much of the early knowledge of these effects 
was obtained by observing only the positron energies. In later 
years both spectrometers gained the capability of simultane- 



ously observing the positron and electron pairs. Both spec- 
trometers were constructed assuming that the entity emitting 
the electron positron pairs had its initial momentum along the 
beam axis. In this case the pairs should have been observed 
in coincidence. However, it was observed, especially by the 
EPOS group that there could be quite significant delays in the 
arrival of both the positron and electron. It is now clear that 
a significant cause of this could be the fact that the scattered 
ions could be not moving along the original beam direction 
when the free electron positron pairs were emitted. 

Table II gives the fit of theory to experiment for the positron 
energies only. For this comparison the application of the the- 
ory is very simple, the positron energy is assumed to be half of 
the energy of the pair The comparisons are attempted only for 
the major peaks observed by experiment. Each line of the ta- 
ble includes: (1) the experiment identifier, (2) the experimen- 
tal peak value, (3) the theory peak value for 9^ = 7r/4, (4) the 
sign (±) for the je solution, (5) the transition A : X ^ Y, and 
(6) the value of 0^ which gives the experimental peak value. 
IN some entries n Table II, where the energies are very close, 
two cases are listed in one row. For example, in the first en- 
try [/ + Cm has two transitions U:K— J^Ll' andU:K— >L2' and 
these are combined to read as U:K— >Lr,L2' and the corre- 
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spending energies and opening half angles are listed with a 
comma. For some of the higher energy transitions the energy 
levels would have to be higher than K^ LI, L2 and are referred 
to as Z— >Z' . In these situations the values are not very differ- 
ent and the transitions are listed as, for example, U, Th:Z— >-Z' 
to indicate that either ion could be involved. 



V. SUMMARY 

The D2 alternative, the "path not taken" by Dirac and oth- 
ers, appears to correct a number of inconsistencies that ap- 
peared as a result of the Dl alternative. In the D2 alterna- 
tive, the reflection symmetry about zero energy is preserved 
in the presence of potentials. The Dirac equation including 
interactions with a vector potential is shown to have charge 
conjugation invariance. Appropriate Ehrenfest theorems are 
recovered. The minimal coupling substitution for the D2 al- 
ternative is found to be 



pt' ^pt^ ^ sgn{E)eA''{x, t), 



(108) 



here as in the rest of the paper, SI units are being used. Some 
surprising results were found for constant potentials in the D2 
alternative. Negative potentials for positive energy states can- 
not be deeper than — m. Positive energy and negative energy 
levels cannot cross the zero energy axis. 

Positive energy potentials Vq that could make up barriers 
do not have any upper limit on Vo . The allowed energy for a 
wave function is only for E > Vq + m for positive energies 
and E < — (Vo + m) for negative energy states. There are 
no propagating wave states allowed in the interval (— Vq, Vq). 
The absence of states below the barrier height Vq completely 
suppresses the Klein tunneling paradox. In the D2 alternative 
the physical interpretation that high static barriers or potential 
wells induce particle hole production is no longer necessary. 

There is no Klein paradox or tunneling in the D2 alterna- 
tive. As already indicated by Dragomanf21], the usual in- 
terpretation that the Klein paradox for a static potential indi- 
cates the production of electron positron pairs is not needed 
theoretically in Dl. Neither is it needed in D2. The more 
troubling Klein tunneling is excluded in alternative D2 and 
so there is little theoretical support for the interpretation of 
the Klein paradox as evidence for pair production by static 
potentials. Pair production and annihilation due to dynami- 
cal effects obviously remains an important part of relativistic 
dynamics. 

In D2 the fact that the charge current density is propor- 
tional to the momentum completely replaces the Dl require- 
ment that negative energy states were necessarily part of any 
wave function, that confinement implies greater negative en- 
ergy components, and that the current density exhibits zitter- 
bewegung. All of these Dl properties are logically unnec- 
essary in D2 and thus unnecessary as part of the relativistic 
theory. 

There is strong evidence that the observations of free e~e+ 
pairs in heavy ion scattering experiments can be rationalized 
in the D2 alternative. 



The possibility of pulling positive energy states down into 
the vacuum is forbidden by the symmetry induced by the 
charge conjugation invariance of the D2 alternative. This is, at 
least, consistent with the experimental failure to find unstable 
vacua in the heavy ion scattering program. 

Finally, there are probably many more questions that re- 
main to be explored about the D2 alternative. The discussion 
here has been very preliminary, but promising. Our prelimi- 
nary stance is that it appears that Dirac and the physics com- 
munity should have taken the D2 alternative in the first place. 
And that could have made all the difference lUTII . 
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Appendix A: Derivation of the Lorentz Force Law 

In this appendix the equation of motion for the generalized 
momentum is calculated for the Dirac Hamiltonian in Alter- 
native Dl. This derivation follows closely the treatment by 
MerzbacherflSl except in these equations the charge e is the 
actual charge of the electron, that is, e < 0. The Hamiltonian 
is 

H = ca- {p~ eA) - e<P + p-mc'' . {A - 1) 

Evaluating the commutators in the equation of motion yields 



— (p - eA) = +e{E + vx B), 
at 



(A -2) 



as expected. In order to re- write this for a more direct compar- 
ison with classical expressions the following identity is used 

a{H -e<P) + (H-e<P)a = 2c{p-eA), {A~3) 

then the equation becomes 

— (— (t;(iJ + e<P) + {H + e^)v)) = -e{E + vxB). 

(A-4) 
If the expectation values of this expression are taken to 
recover the non-relativistic Ehrenfest classical equations of 
motion, the expectation values must be taken with respect to 
some states of energy either positive or negative. In this case 
we need to substitute for the following 



{H - e<i>) « sgn{E)mc^, 



(A -5) 



where E is the energy of the state, into the equation of motion, 
yielding 



sgn{E) — {{mv)) = e{E + v x B). 



(A -6) 
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This conclusion within Dl yields a classical expression for 
the momentum, but the Lorentz force law direction depends 
on the sign of the energy of the state determining the expecta- 
tion value. 

Within the alternative D2 the discussion begins with the D2 
minimal coupling substitution 



pt^ ^pf^ - sgn{E)eAi' 



{A -7) 



which will generate a new factor of sgn{E) on the right hand 
side which will cancel out the same factor on the left hand 
side, recovering the appropriate Ehrenfest theorem for the 
classical Umit 



dt 



{{mv)) = e{E + v x B). 



(A -8) 



This is exactly the Lorentz equation for particle with charge e. 



Appendix B: Derivation of tlie Larmor Spin Dynamics 

The role of the spin in the alternative Dl is demonstrated 
by the equation of motion of the Matrix operator S where 



a 
a 



{B-l) 



and a are the Pauli Spin matrices. This discussion again fol- 
lows closely from Merzbacher[18|. The equation of motion 
is now for a Hamiltonian with only a vector potential 

H = ca-{p-eA) + Pmc^ {B - 2) 

and proceeds from 

In preparation for taking expectation values with Hamiltonian 
eigenstates, consider 



H—- + -—H = -2ecS X B. 



dt 



dt 



(5-4) 



In the lowest order expectation value we can substitute H = 
sgn{E)mc^ and this yields for the expectation values 



dt mc 



{B-5) 



which again has the dependence on the sign of energy of the 
eigenstates used in the expectation values. 

In the alternative D2 the magnetic field will end up with a 
factor of the sgn{E) on the right hand side and this will re- 
cover the classical Larmor equation for the spin in a magnetic 
field 

^ = -^(E)xi?. iB^6) 

dt mc 



Appendix C: Metastable State Decay Effect on Experiments 

The decay of these metastable excited ions after the colli- 
sion with the target can be studied approximately by evalu- 
ating Feynman diagrams. There are at least three processes 
that are likely to be most important. The first of these is the 
rate i?^ of spontaneous emission of a photon as the occupied 
positive energy electron drops into the empty negative energy 
state. The second rate R^p is the absorption of an initial free 
e^e+ pair and the decay of the bound e~e~^ state into a free 
electron -positron pair This free e~e~^ pair is what was de- 
tected in most of the heavy ion scattering experiments. A 
third rate Rind would be the induced decay of these bound 
e~e~^ pairs because of other ions, pairs, photons, and time 
varying electric fields in the scattering region after the target. 
This induced decay should be proportional in lowest approxi- 
mation to the square of the current since it reflects interactions 
between beam ions and their associated post scattering con- 
stituents. The other two rates should be proportional to the 
beam current. This means that the production of free e^e+ 
pairs should be proportional to 



R<p + Re 



Ri: 



Since every term in this ratio is at least proportional to the 
beam current, we can factor out one power of this current and 
get an expression for the production cross section for e^e+ 
pairs that has the following structure. 



cr^ 



11 J 



(C-1) 



where 77 J = Rind/ Rep is a parameter expressing the ra- 
tio of the induced emission term divided by the spontaneous 
emission term, J is the beam current, and Xq = {R4, + 
Rep) /Rep and is independent of the beam current. 

Let us now examine the effect of this expression on the time 
it takes to collect a statistically significant signal in a scatter- 
ing experiment. In carrying out this analysis, let us express 
the cross section for the background e^e+ pairs as 0-2 ■ In 
order to establish a statistically significant signal we need the 
ratio of the signal over the standard deviation to be some pre- 
determined number, z , which can be expressed as 



(T^pJL 
\fo^2Jt 



(C-2) 



where t is the time the counting experiment must run to 
get z . In the usual scattering experiment (essentially the as- 
sumption behind the APEX papers), there is no factor reflect- 
ing competing channels for the supply of the metastable states 
that give rise to the e^ e+ pairs, and we can replace u^p by u\ 
and get the following expression for the time t that it takes to 
achieve statistical significance. 



^z'^a2/{(Jlp 



to 



(C-3) 
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where for the following comparison we have expressed the 
formula in dimensionless variables x ~ rjj and r = t/to. If 
we now write in these units the standard expectation for an 
experiment, we get 



(C-4) 



which reflects the experimental wisdom that increasing the 
beam current decreases the counting time to achieve statis- 
tically significant results. 

However, the same expression for the case where 
metastable bound e^e+ states are competing with the spon- 
taneous and induced emission of photons yields the following 
form for statistically significant observation times : 



{Xo 



iC-b) 



As can be seen by the functional dependence on x this func- 
tion has a minimum at Xo and then increases as the beam cur- 
rent (x) increases. In this case, large currents will actually 
increase the counting time, and there is an optimal range of 
beam currents which gives the shortest counting time for ob- 
serving the emitted free electron positron pairs from the decay 
of the metastable states. 

Two very different behaviors are to be expected for the 
counting time to significance. If there is no metastable states 
that can be induced to decay, increasing the beam current, al- 
ways decreases the counting time. However, in the presence 
of a metastable, inducible, intermediate state, there is an opti- 
mal counting time and increasing the beam current increases 
the counting time considerably and can appear to overwhelm 
the peaks if not enough time is used for counting. 
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